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Transport properties of an interacting triple quantum dot system coupled to three leads in a tri- 
angular geometry has been studied in the Kondo regime. Applying mean-field finite-U slave boson 
and embedded cluster approximations to the calculation of transport properties unveils a set of rich 
features associated to the high symmetry of this system. Results using both calculation techniques 
yield excellent overall agreement and provide additional insights into the physical behavior of this 
interesting geometry. In the case when just two current leads are connected to the three-dot system, 
interference effects between degenerate molecular orbitals are found to strongly affect the overall 
conductance. An 5 = 1 Kondo effect is also shown to appear for the perfect equilateral triangle sym- 
metry. The introduction of a third current lead results in an 'amplitude leakage' phenomenon, akin 
to that appearing in beam splitters, which alters the interference effects and the overall conductance 
through the system. 

PACS numbers: 73.63.Kv, 72.10.Fk, 72.15.Qm 



I. INTRODUCTION 

One of the most important and exciting aspects of 
molecular physics nowadays is the study of electronic 
transport properties of natural and/or fabricated struc- 
tures at the nanoscopic scale. In nature, molecules 
can couple to the external environment (electron reser- 
voirs) through extended orbitals, which permit conduc- 
tion electrons to hop in and out of the molecule. Fabri- 
cated molecules can be made by coupled quantum dots 
(QDs) with discrete energy levelsji'^'^ Depending upon 
the strength of the coupling between the QDs, they can 
behave as a molecule with extended orbitals, which can 
further couple to external electron reservoirs. The con- 
finement of electrons inside this artificial 'molecule' pro- 
duces strong Coulomb interactions)^ which may give rise 
under suitable conditions to Kondo physics for temper- 
atures below a characteristic crossover temperature, the 
Kondo temperature, Tk^'^^ In the simplest picture of 
this regime, for T <^Tk, the system forms a singlet state, 
created by the screening of the localized spin by the con- 
duction electrons in the external reservoir. Since its first 
observation in QDs in 1998,^ the attention generated by 
the Kondo effect in these structures has led to an explo- 
sion in experiments and theory. For example, multiple 
QD systems have become platforms for the theoretical 
and experimental development of sophisticated arrange- 
ments in order to access the rich phenomenology of the 
Kondo problem, including non-Fermi-liquid behavior and 
quantum critical pointsi^ 

In this paper, we study the transport properties of 
a triple quantum dot (TQD) system in the 'molecular 
regime' [with strong interdot couplings; see Fig. [TJa)] 
in two distinct situations: Firstly, just two QDs are at- 



tached to independent electron reservoirs. Secondly, each 
QD is connected to an independent electron reservoir. In 
the latter case, we focus our attention on the conduc- 
tance of the system through two of the three terminals 
(the same two used to measure conductance in the first 
case) . We are particularly interested in understanding in- 
terference effects, especially the role played by the third 
lead in the propagation of electrons along the different 
trajectories. 

Despite significant advances in the understand- 
ing of Kondo physics in doubl o^°i^^'^^i^^i^"^ and 
tripl o^^i^^i^'^i^^i^^i^" QD structures made in the last few 
years, there are still important aspects of the problem 
which deserve to be studied in detail. For example, based 
on a suggestion by Zarand et al.,^^ one may ask if an 
SU(4) Kondo regime may be experimentally attained in a 
TQD geometry. In addition, the unprecedented control of 
parameters in these multi-dot structures opens the pos- 
sibility of observing quantum critical points and their as- 
sociated non-Fermi-liquid ground states. Although many 
of these have been theoretically identified, the very de- 
manding experimental constraints required have resulted 
in only a few successful experimental realizations.^ Fur- 
ther motivation to study TQD systems comes from the 
proposal by Saraga and Loss that these structures could 
be used to produce spatially separated currents of spin- 
entangled electronsi^ Experimentally, however, only few 
groups have reported work in these systemsi^^iSliS Most 
of these studies have been in the Coulomb blockade 
regime, and one of the works reports that a TQD de- 
vice can act as a molecular rectifier. ^■^ 

Zitko and Boncaii^ have recently studied theoretically 
a TQD system connected in series to two leads. They 
have found that for a certain range of inter-dot hopping 
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parameters the system crosses over from a Fermi-liquid 
to a non-Fermi-liquid regime in a wide interval of tem- 
peratures. In a subsequent paper,— using the numer- 
ical renormalization group (NRG), these authors ana- 
lyze a large number of phases for a system similar to 
the one we discuss in this paper. Notice, nevertheless, 
that there are important differences between their sys- 
tem and ours: The majority of the phases analyzed in 
detail in Ref. use a Kondo Hamiltonian for the dots. 
In that case, a crossover was predicted between the two- 
impurity and the two-channel Kondo-model non-Fermi- 
liquid fixed points. Their analysis of Anderson- impurity 
QDs brings their work closer to ours. However, they re- 
strict their study mainly to regimes where the inter-dot 
hoppings are considerably smaller than the coupling to 
the leads, which is exactly the opposite regime we treat 
in our work. In addition, they only analyze results close 
to half-filling, while we consider all fillings. Finally, and 
more importantly, they did not analyze the very impor- 
tant influence of the third contact, which is one of the 
important results of the work we report here. 

The physics of this arrangement of QDs has also been 
the subject of other theoretical works . ^''"i^^ '^'^'^'^'^ In 
particular, a situation where the system may present in- 
teresting, but rather complicated behavior, is in the fully 
symmetric case, i.e., when all inter-dot hoppings are the 
same and each QD is equally connected to an indepen- 
dent conducting band (in that case, the system has equi- 
lateral triangle symmetry). It is reasonable then, if the 
inter-dot hoppings are much smaller than the intra-dot 
Coulomb repulsion, to expect the TQD system to present 
a spin frustrated regime, as anti-ferromagnetic arrange- 
ment between electrons sitting in different QDs is not 
possible. Indeed, through the use of conformal field the- 
ory and NRG calculations, Ingersent et al. were able 
to characterize a novel, stable, frustration-induced non- 
Fermi-liquid phase for a three-impurity Kondo modeljSS 
It should be noted that this is not the regime treated 
in the current work. Here, we concentrate in the regime 
where the inter-dot hoppings are of the same order of 
magnitude as the intra-dot Coulomb repulsion, and al- 
ways larger than the coupling to the leads (the 'molecular 
regime'). 

Considering this rich theoretical context, it is impor- 
tant for our objectives to be clearly spelled out. They 
are three-fold: First, since it is important from an ex- 
perimental point of view to analyze the charge fluctua- 
tions in the QDs as a function of the gate potential, and 
as most of the previous work mentioned above uses the 
Kondo model to represent the QDs, we will model the 
system using the Anderson impurity model to describe 
each quantum dot. Second, we carefully analyze the con- 
ductance vs. gate potential results in a regime where the 
inter-dot couplings are larger than the coupling to the 
leads, i.e., in the molecular vegivoB. Although this regime 
excludes other interesting phases in this system analyzed 
before we believe that the molecular regime can be 
experimentally more accessible and therefore very rele- 



vant. Third, we analyze in detail the effects created by 
the introduction of a third electron reservoir, which is 
connected to the 'free' QD, i.e., the QD which is not 
connected to either of the reservoirs used to measure the 
conductance. 

We study this system by calculating the appropriate 
propagators to obtain the charge, the local density of 
states at the dots, and the conductance, using two differ- 
ent approaches: a finite-U slave boson formalism devel- 
oped in the mean- field approximation (FUSBMF)^ and 
the Embedded Cluster Approximation (EC A) where 
one diagonalizes a small cluster containing the dots, and 
then embeds it into the leads through a Dyson equa- 
tion. These two completely different approaches provide 
a similar description of the physics of the TQD struc- 
ture. Note that some of the results shown here were 
obtained using a recent variant of the EGA method, the 
Logarithmic-Discretization Embedded Cluster Approxi- 
mation (LDEGA).^''^ In this variant, the non-interacting 
electron band is discretized logarithmically (a la NRG), 
which leads to much faster convergence with cluster size. 
The band discretization provided by LDECA is neces- 
sary in two circumstances: i) When finite-size effects pre- 
clude EGA from converging to the correct ground state, 
or when that convergence is too slow; and ii) When one 
wants to calculate quantitatively accurate local density 
of states (LDOS). We will clearly indicate when either 
method (EGA or LDECA) is used. 

This paper is organized as follows: In section II, we 
specify the model used to represent the TQD and we 
briefly describe the numerical methods used (FUSBMF, 
EGA, and LDECA). The results obtained for the case 
where the TQD is coupled to two reservoirs are discussed 
in section III. The change in the transport properties 
caused by the introduction of a third lead attached to 
the TQD is discussed in section IV. Finally, in section V, 
we present the conclusions. 



II. TQD MODEL AND NUMERICAL METHODS 

The TQD system studied in this work is schematically 
represented in Fig. [IJa). The full Hamiltonian can be 
written as 

H = HqDs + Hieads + Hqus-leads, (1) 

where Hqds describes the isolated TQD system, Hieads 
the two (or three) independent leads and Hpos-ieads es- 
tablishes the contacts between the dots and the leads. 
Explicitly, we have 

i = A,B,C i = A,B,C 

+ H.C.] , (2) 
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FIG. 1: (color online) (a) Schematic representation of triple 
quantum dot system coupled to leads. As indicated by the 
arrows, we will calculate the conductance between the left (L) 
and right (R) contacts, (b) Energy diagram (in units of t^) 
for the molecular orbitals [as defined in Eqs. (5)] as function 

of ti/t-j. 



where dj^ (dia) creates (annihilates) an electron with spin 
i*'* QD and energy controlled by the gate poten- 
= d\^dia is the occupation number operator, 
and U is the Coulomb repulsion energy for double oc- 
cupancy in a QD. The leads, modeled as semi-infinite 
chains, are represented by the Hamiltonian, 
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tial Vgi, n. 
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H.c. 



(3) 



where cj ^ {ci a-) creates (annihilates) an electron with 



spin a in the i site of the 7 lead, and t is the kinetic 
hopping between first neighbor sites. Finally, the con- 
tacts between the QDs and the leads are established by 
the Hamiltonian 



QDs-leads - '^\tl{d\„Ci^c, + d^u^Ci^^a) 

+ ^2dLcip^ + H.c.l , (4) 



where ci^o- (7 — L, R, P) annihilates an electron in the 
first site of the 7*'' lead. Note also that for all cases 
studied here, a left <-> right symmetry applies. 

In order to gain some intuitive understanding of the 
TQD system, let us consider the non-interacting case 
first. Let us further assume that ^1=^2 = 0, so that 
the QDs are completely disconnected from the current 
leads. This is the 'atomic limit' of the model and can be 
solved exactly. In that case, the Hamiltonian eigenvalues 
are {Vgi ^Vg), 
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E2 
E3 
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Vg-U 



8tl 



(5a) 
(5b) 
(5c) 



The un-normalized states are given by ([A,B,C]), 
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[-1,0,1] 
[1,-— - 
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^3 



<1 + 8*1,1] 
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2*3 



7^\Al + 8ii,l]. 



(6a) 
(6b) 
(6c) 



Borrowing the terminology from molecular physics, or- 
bitals I ipi), I 1P2), and I -03) will be denoted as bonding, 
non-bonding, and anti-bonding, from now on. For the 
particular case of t4 = t3 = t" , Ei = E2 = Vg — t" and 
E3 = Vg + 2t", the system has a doubly-degenerate state. 
In this case, the eigenvalues Ei, E2, and E3 correspond 
respectively to the orbitals 



^1) 
^3) 



1 

1 

71 



1,-2,11, 



-1,0,1], 



(7a) 
(7b) 
(7c) 



The degeneracy results from the symmetry of the sys- 
tem. In group theory language, it is associated to a two- 
dimensional irreducible representation of the Ca^ sym- 
metry group. Note that, obviously, each orbital is also 
SU(2) symmetric, therefore, at zero field they are doubly 
degenerate regarding the spin orientation. For the full in- 
teracting Hamiltonian, these orbitals hybridize with the 
conduction electron band, renormalizing the eigenvalues. 
Although this is a simplified picture, it helps to under- 
stand the transport properties of the system. We will 
show that the degenerate orbitals have an important in- 
fluence in the conductance of the interacting case. 

In the interacting case, we are mainly interested in the 
Kondo regime and will study in detail how the symme- 
try (with its associated degeneracy) affects the transport 
properties. As mentioned above, the system is analyzed 
applying the FUSBMF, EGA, and LDECA methods. All 
three methods allow the calculation of the Green's func- 
tions. We can easily calculate the charge at each dot and 
the total conductance of the system, which are respec- 
tively given by 
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lm{GUto)}f{u;)d^ 



and 



Qrp ^ t-^^pR{€F)pL{<^F) I Glr{^f) 



(8) 



(9) 



where Gfj is the local Green's function of the QDs, f{iS) 
is the Fermi function, Pl(r){^f) is the density of states of 
the left (right) lead's first site and GLRiep) is the Green's 
function that propagates an electron from the left to the 
right lead, all calculated at the Fermi energy e^- The 
expression for the conductance can be derived from the 
Keldysh formalism'^— and is equivalent to the Landauer- 
Biittiker formula for the non-interacting case. 
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A. Finite-U slave bosons mean-field approximation 

In the FUSBMF approach, one enlarges the Hilbert 
space by introducing a set of slave boson operators Ci, 
Pia and di {i = A,B,C), and replace the creation (djg.) 
and annihilation {dta) operators in the Hamiltonian by 
d\^z\^ and Zi^dia, respectively. Following Kotliar and 
Rukenstein, the operator z takes the forro^ 

Z,^ = [1 - d\di - p\„P^crY''^[e\'P^a + p\^di] 

A^~e\h~plp^-.Y'^■ (10) 

Notice that the bosonic operators di and ei do not carry 
spin index. The enlarged Hilbert space is then restricted 
to the physically meaningful subspace by imposing the 
constraints 

P^ = e\e, + ^ p\j^, + dH -1 = (11) 

a 

and 

Qza = dl^dia - pl^Pia - = 0. (12) 

These constraints are enforced by introducing them into 
the Hamiltonian through Lagrange multipliers A^^^ and 

XlJ . The constraints pT|) force the dots to have empty, 
single or double occupancy only, and ^2]) relates the bo- 
son to the fermion occupancy. In the mean-field approx- 
imation, the boson operators Ci, pia and di (and the cor- 
responding Hermitian conjugates) are replaced by their 
thermodynamical expectation values = (e^) = (e|), 

Pia = (Pia) = (pla) and di = (di) = (dj). These expec- 
tation values, plus the Lagrange multipliers, constitute 
a set of parameters to be determined by minimizing the 
total energy (H). In principle, a set of seven selfcon- 
sistent parameters are needed for each dot. Although 
our system has three dots (which would require a total 
of 21 parameters), we take advantage of the symmetry 
of the configuration, since QDs A and C are symmetric 
with respect to B. Note that this symmetry imposes no 
additional constraint on the parameter values. In con- 
trast with a previous implementation of this methodjii 
our approach allows for a more complete and versatile 
description of the system in terms of its structural pa- 
rameters. In particular, it can describe the transition 
from non-symmetrical to highly symmetrical regimes as 
the interdot parameters are changed. 

In the mean-field approximation, we can obtain self- 
consistent expressions for the Green's functions. Then, 
using Eq. ([5]) and ^ , we can calculate the charge at each 
dot and the conductance. We note that all the calcula- 
tions with FUSBMF are implemented on the individual 
QD basis, while EGA and LDEGA utilize the molecular 
basis. 



B. Embedded Cluster Approximation 

The EGA method^i^ relies on the numerical determi- 
nation of the ground-state of a cluster with open bound- 
ary conditions. In the following, we briefly sketch details 
of the method. 

The EGA method tackles the impurity problem in 
three steps. First, the infinite system is naturally cut 
into two parts: one part C (the cluster) contains the 
interacting region plus as many noninteracting sites of 
the leads as possible, and a second part R (the 'rest'), 
consisting of semi-infinite chains positioned at left and 
right in relation to the cluster C. The number of sites 
in C is denoted by A'ed- Second, Green's functions for 
both parts are computed independently: current imple- 
mentations of EGA utilize the Lanczos method^ to cal- 
culate the interacting Green's function of the interacting 
region, while those of the part R, being noninteracting, 
can be computed exactly as well. In a final step, the ar- 
tificially disconnected parts are reconnected by means of 
a Dyson equation, which dresses the interacting region's 
Green's function. This step, the actual embedding, is 
crucial for capturing the many-body physics associated 
with the Kondo effect. Moreover, although the clusters 
that can be solved exactly by means of a Lanczos rou- 
tine are rather small, being of the order of TVed ~ 12 
sites only, the embedding step successfully compensates 
for that by dressing the cluster Green's function and ef- 
fectively extending the many-body correlations, induced 
by the presence of the impurity, into the semi-infinite 
chains R. Obviously, strongly correlated regimes which 
depend on extremely low energy scales will be difficult to 
treat with the embedding procedure, although, as men- 
tioned above, great progress has been made lately in this 
respect by introducing a logarithmic discretization pro- 
cedure into the algorithm. 

We now provide further detail on these steps. The 
Hamiltonians of the left and right semi-infinite, tight- 
binding chains, i.e., the noninteracting R part, are de- 
scribed by 

— oo 

^^sc-L = ^ {c\a'^l-la + H.c); 

oc 

Hsc-K = -t J2 + H.c), (13) 

l=NEr> + l,a- 

where in this notation, the sites labeled by i = 1, . . . , Neb 
are inside the cluster C. The semi- infinite chains are con- 
nected to the cluster by the following term: 

i/hy = -V[c\„Coa + C^N^n'y^^En + l'y] + H.C, (14) 

where V — t is the hopping in the broken link, connect- 
ing parts R and C, used in the embedding procedure. 
The Green's function for the cluster C and for the semi- 
infinite chains are calculated at zero temperature. Fixing 
the number of particles m and the z-axis projection of 
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the total spin, S^^^^^, the ground state and the one-body 
propagators between all the clusters' sites are calculated. 
For example, 5^™''^*°*°''', the undressed Green's function 
for the cluster, propagates a particle between sites i and 
j inside the cluster. For the noninteracting, semi-infinite 
chains, the Green's functions gQ and g§^^^i at the sites 
and -/Ved + 1, located at the extreme ends of the semi- 
infinite chains, at left and right to the cluster, can be 
easily calculated as well. 

The Dyson equation to calculate the dressed Green's 
function matrix elements 6',-'"''^'°'''''' can therefore be 
written as 



G 



('"'■S'total) _ „(™'-Stotal) I A"^-'Stat^l) 



2J 



+ 9i,i 



-'total'' ^ ^V"^''-'total'' 



r'('"''5'?otai) 

^(".•S^otal) 



9Nev + 1 ^ ^Neu,] 



(15) 

(16a) 
(16b) 



where V ^ as mentioned above, is defined according to 
i/hy Eqs. US]) and (|16p correspond to a chain approx- 
imation, where a locator-propagator diagrammatic ex- 
pansion is usedj^i^ Note that EGA is exact in the case 
of C/ = 0. 

As mentioned before, the calculation of the propaga- 
tor 5^™'^'°'°''' requires that fixed quantum numbers m 
and iS'^q^jji be used. However, after the embedding pro- 
cedure, these quantum numbers are not good quantum 
numbers for the cluster anymore. Therefore, we have to 
incorporate processes into the EGA method that allow 
for charge fluctuations in the cluster C. To accommo- 
date this requirement, different implementations of EGA 
have been devised, either by including different spin mix- 
ing strategies^i^ or by moving the Fermi energy of the 
leads. ''I 

The spin mixing proceeds as follows. First, a cluster 
Green's function with mixed charge is defined through 

{m+p.pSL^,,) /-, N (m,0) , (™+l>'Stotal) 

9ij (1 -P) 5ij ' + P 9lj , (17) 

where p takes values between and 1, and we are as- 
suming that m is even, in which case, the corresponding 
•^totai = 0- In addition, note that for the cluster with 
'total takes values ±1/2. The matrix ele- 



charge m + 1, S", 

ment (;|™+P'P'^tot'>i^ corresponds to a situation where the 
charge in the cluster is between m and m -f 1. The to- 
tal charge in the cluster, before embedding, can be easily 
calculated as 

qPStat^i (p) = (1 - p)m + p{m + 1) ^ m + p. (18) 

Using Eqs. (fT5)) and (|16p. the dressed Green's function 
^(m+p,pSt„t^i) obtained, and from this result, the charge 
in the cluster can be calculated: 

Im { 



IT 



""'(u;)}dw, 
(19) 



where Ep is the Fermi level. The value of p is calculated 
self-consistently, satisfying 



QPStotal(p) = qPS^=tal(p). 



(20) 



If there is spin reversal symmetry, e.g., no magnetic field 
is applied, one can calculate the total Green's function 

as 



E 



G 



{m+p.pS', 



(21) 



=±1/2 



where p satisfies Eq. (j20p . It is important to empha- 
size that the charge fluctuations taken into account by 
Eq. p7)) are the ones between the cluster and the rest 
of the system, and not just the ones at the interacting 
region described by i?int- The latter ones involve a very 
localized neighborhood of the dot and as a consequence, 
are typically already well described on isolated clusters 
only. Finally, it is noteworthy to point out that the self- 
consistent solution for the charge mixing parameter p 
is either or 1 in the Kondo regime, and in particular 
at the particle-hole symmetric point Vg = —U/2 (when 
analyzing a single- QD problem). Therefore, deep into 
the Kondo regime (i.e., at the particle-hole symmetric 
point), no charge mixing takes place at all, and very lit- 
tle charge mixing occurs in a window of gate potential 
around the particle- hole symmetric point. The parame- 
ter p will start to take a finite value (note that < p < 1) 
as the gate potential drives the system into the mixed- 
valence regime. The purpose of the charge mixing is thus 
mainly to smooth out the transition from an N electron 
to an ± 1 electron ground state, which for the bare 
cluster is a crossover between ground states with differ- 
ent number of particles. 

As mentioned above, some of the calculations were 
done using the LDEGA method, which is an impor- 
tant extension of EGA. In it, to obtain a better descrip- 
tion of the low energy physics of the system, the non- 
interacting band is logarithmically discretized. All the 
procedure described above remains the same, but the 
band discretization allows a much faster convergence to 
the Kondo regime with cluster size. A full description of 
LDEGA can be found in Ref. [M 



C. Numerical results 

In order to study the conductance as a function of the 
parameters of the system, we use the leads' hopping cou- 
pling t as the energy unit {t = 1), and set the Fermi 
energy to zero {ep — 0). All the results will be shown 
for zero-temperature. In the strong interdot coupling 
regime, i.e., ts^t^ > ti,t2, the individual QDs levels mix 
into three molecular orbitals which are coupled to the 
leads. The energy of these orbitals can be controlled by 
gating the local energy states and by varying the hopping 
matrix elements ratio ^4/^3 [see Fig. [IJb)]. Their widths 
depend upon the coupling to the conduction bands, i.e., 
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FIG. 2: (color online) Conductance as function of gate po- 
tential, Vg = VgA = VgB = Vgc , foi varlous values of t^, as 
indicated in the panels. (Blue) solid dots indicate LDECA 
results for clusters with 11 sites; (red) 'x' signs are FUSBMF 
results. Notice in panel (a), where ti = 0, that the number 
of electrons in the TQD system is indicated for each one of 
the valleys, where the Kondo effect is absent. The structure 
evolves from three to two peaks as the system goes from three 
dots in series [f4 = 0.0, panel (a)], to isosceles \ti < ta = 0.5, 
panels (b) and (c)], and finally to equilateral triangle symme- 
try [t4 — ts, panel (d)]. For all panels, U = 1.0, ti = 0.45, 
t2 — 0.0, and ts — 0.5. The very small discrepancy between 
LDECA and FUSBMF in panel (d) comes from a finite-size 
effect in the LDECA results. 



ti and t2- In order to study the contribution of each in- 
dividual orbital to the conductance, we make them suf- 
ficiently far apart from each other. To do so, we take 
^3 > ti,t2- In particular, in section III, we set ti — 0.45, 
t2 — 0.0, = 0.5, and U = 1.0 (which are the same 
parameters used in Ref. ITtI ). and in section IV we set 
h = 0.2, 0.0 < t2 < 0.2, ^3 = 0.4, and U = 0.5^^ In 
section III, we vary to manipulate the symmetry of 
the system (0.0 < t4/i3 < 1.0). In section IV, besides 
the same variation of ti as in section III, we also ana- 
lyze what is the effect of varying t2 (0.0 < t2/ti < 1.0), 
i.e., we verify what is the effect of adding a third lead 
[connected to QD B, see Fig. [TJa)] to the TQD system. 



III. TQD CONNECTED TO TWO LEADS {t2 = 0) 

A. TQD in series (ti = 0) 

Initially, taking advantage of the flexibility of the nu- 
merical methods used, we analyze the conductance when 
the three QDs are aligned in series [ti = 0) and cou- 
pled to two leads only (i.e., we make t2 = 0). Fig. [2] 
shows the conductance as a function of the gate potential 



Vg = VgA = VgB = VgQ ■ pauel (a), for t^ = [(red) 
'x ' signs indicate FUSBMF results, and (blue) solid dots 
display LDECA results] , the three molecular orbitals are 
equally separated by an energy value proportional to t^. 
As will be shown below, this originates from three Kondo 
peaks (occurring at different Vg values) associated to each 
molecular orbital. 

The t4 — case corresponds to the molecular regime 
reported in Ref. fisl (with, as mentioned above, ia > ti). 
For t4 7^ 0, the system transforms into a triangular con- 
figuration, which can be compared to the system studied 
in Ref. [l3. In the first case {t^ = 0, panel (a) in Fig. [2]), 
the three peaks in the conductance occur at gate poten- 
tial values where there is a change in the occupation of 
the different molecular orbitals. When the bonding or- 
bital hosts one electron (for Vg/U — 0.5), the system is in 
the Kondo regime and the characteristic Abrikosov-Suhl 
resonance of this regime creates a path for the electrons 
to cross from the left (L) to the right (R) lead [this will be 
more clearly demonstrated in Fig. [3]^d)]. Decreasing Vg, 
we find a valley corresponding to the accommodation of 
a second electron in the bonding state, creating a singlet 
that destroys the Kondo effect. The middle peak corre- 
sponds to the presence of a third electron in the system, 
now siting in the non-bonding state, since the bonding 
state is full. Again, the conductance peak reflects the 
Kondo resonance at the Fermi level due to an unpaired 
electron that is anti-ferromagnetically correlated with the 
conduction electrons (see Fig. [3)Jb), and discussion be- 
low). The following valley and the third peak are a con- 
sequence of the suppression of the Kondo effect due to 
double occupation of the non-bonding orbital and the un- 
paired electron in the anti-bonding orbital, respectively. 
Finally, the final drop in the conductance results from 
the destruction of the Kondo effect due to the sixth elec- 
tron entering into the system. The electron occupancies 
at the conductance valleys are indicated in panel (a). 



B. LDECA LDOS results for molecular orbitals 

it4 = 0) 

Before analyzing the effect of introducing a finite t^, 
we want to show LDOS results at the particle-hole sym- 
metric point (Vg = —U/2) [Figs. [3{a) and (b)], and at 
Vg = U/2 [Figs. El^c) and (d)] for the t4 = 0.0 curve in 
Fig. IHa). In the upper left panel in Fig. [3l we have the 
LDOS for each QD for Vg = -U/2. As expected, QDs A 
and C have the same LDOS (dark (black) solid curve), 
and the peak at w = is indicative that they participate 
in a Kondo effect. Indeed, a clear Kondo peak can be 
seen at the Fermi energy Ep (uj ~ 0.0) for QDs A and 
C, while the LDOS for QD B has a gap at Ep [dashed 
(red) curve in panel (a)]. This will be important later 
on to understand the results when the upper lead (P) is 
coupled to QD B (for finite tz). The LDOS for the ap- 
propriate orbital states for this configuration (ti = 0.0, 
and Vg = —U/2) is shown in the lower left panel. The 
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co/U (O/U 

FIG. 3: (color online) LDOS results using LDECA (9 sites) 
for ti = 0.45, t2 = 0.0, ts = 0.5, U = 0.0, and U = 1.0. (a) 
and (c): LDOS for each QD. (a) and (b) are calculated at 
Vg = -(7/2, while (c) and (d) at Vg = U/2. Note that, in (a), 
QDs A and C are in Kondo and QD B is not, while in (c) all 
dots participate in the Kondo effect, (b) and (d): LDOS for 
each molecular orbital. In (b), just the non-bonding orbital is 
in the Kondo state (dark (black) curve) , as one can easily see 
its Kondo peak at w = 0.0, while the bonding orbital (light 
(green) curve) is almost fully occupied (most of its LDOS is 
below the Fermi energy) and the anti-bonding orbital (dashed 
(red) curve) is almost empty (most of its LDOS is above the 
Fermi energy). In (d), it is now the bonding state that has 
a Kondo resonance, while the non-bonding and anti-bonding 
states are higher in energy. The inset shows a more detailed 
view of the Kondo peak for the bonding orbital, where the 
small LDOS peaks for the e and e + U states are indicated by 
arrows. 



bonding orbital (gray (green) solid curve) has most of 
its LDOS below Ep^ indicating that it is already almost 
fully occupied. In contrast, the anti-bonding orbital has, 
at this particular gate potential, most of its LDOS above 
Ep and is therefore almost completely empty. The non- 
bonding orbital [which is an antisymmetrical combina- 
tion of QDs A and C, only, see Eq. 6(b)] displays a 
Kondo peak at Ep, which is responsible for the unitary 
conductance seen for ^4 = in Fig. [l^a) at Vg = —U/2. 
As to the rightmost peak in Fig. [D^a), notice, as can 
be seen in the upper right panel in Fig. [Sljc), that all 
three QDs participate in the Kondo effect for Vg = U/2. 
In the lower right panel, one sees the LDOS for the or- 
bital states, now indicating that the bonding state has 
a Kondo peak, while the other two orbitals are nearly 
empty.— The inset shows details of the Kondo peak, for 
the bonding orbital, including the shorter peaks associ- 
ated to e and e + U, indicated by vertical arrows. 



C. Finite ti and interference effects {t4 « ts) 

Now, as t4 increases from to 0.5 (panels (b) to (d) in 
Fig. [2]), the bonding and non-bonding energies become 
closer to each other [see Eqs. (5)], and finally, become 
degenerate for t4 — — 0.5 [see Fig. 1(b)], when the 
system possesses an equilateral triangle symmetry. Note 
that the peaks in the conductance curves shown in Fig. [21 
corresponding to the bonding (rightmost peak) and non- 
bonding (central peak) orbitals, merge into each other 
and the conductance decreases as t^. As shown 

above in Fig. [3l the LDOS of the molecular orbitals pro- 
vides a more clear picture of the Kondo effect than the 
LDOS of each QD. For example, when = 0, the right- 
most conductance peak in Fig. [2la) can be directly as- 
sociated to the Kondo peak in Fig. [3jd) (hght (green) 
solid curve). As mentioned above, the molecular orbitals 
provide a natural description of the Kondo effect when 
the intra-dot hoppings arc larger than the coupling to the 
leads. 

In the same manner that the LDOS for each molecular 
orbital provides important insight into the conductance 
through the TQD system, one can define a 'partial' con- 
ductance Gi through each molecular orbital 'i' (i = L2, 3) 
in the following way (full details are given in Ref . 1441) : 

Gi = ^ [t^gitHp{ep)]^\G,R{ep)\^ (22) 

where gi is the Green's function in the first site of the 
left contact, tu is the coupling of the left lead with orbital 
I tpi), and Gift is the dressed Green's function that moves 
an electron from | V'i) (where i = 1,2,3) to the first site 
in the right contact. For ^4 « ^3 and Vg values such that 
mostly molecular orbitals i = 1,2 are involved in the 
transport of charge {i.e., —0.5 < Vg/U < 0.5, in panels 
(c) and (d) in Fig. [2]), the total conductance Gt can be 
approximated by the equation 

Gt ^ G12 = Gi + G2 + 2VG1G2 cos A</)i2, (23) 
where 

defines the phase-difference between a path that goes 
through orbital | tpi) and a path that goes through or- 
bital \ ip2)- In the case where all three orbitals are con- 
tributing, a simple extension of these equations should be 
used, and it gives results exactly equal to the ones shown 
in Fig. [2l We should note that the Gi functions have a 
characteristic 'width' given by the coupling to the leads 
(and the weight of the orbital at the connecting dot), as 
well as a position dependence on V^, as the energy of each 
orbital shifts with respect to the Fermi energy. 

Equation (|23p shows that when there is no energy over- 
lap between | ipi) and | 1P2) (therefore, no overlap be- 
tween Gi and G2), which occurs when the corresponding 
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orbitals are well separated in energy, the last interference 
term is zero, independently of the value of the phase dif- 
ference. This is essentially the case for the conductance 
results for i4 = in Fig. [H^a). In this case, a simple sum 
of the partial conductances Gi through each molecular or- 
bital is very similar (not shown) to the total conductance 
(and more so as the level separation increases for larger 
values of ts/ii). However, as shown next (see Fig. H]), 
once the molecular orbital levels start to overlap, the par- 
tial conductances G; are no longer simply related to the 
total conductance, as the last term in Eq. now plays 
a role, and its effect will obviously depend on the value 
of A(/)i2. Therefore, the calculation of the partial con- 
ductances, and the phase difference of the corresponding 
Green's functions, provides us with information about 
possible interference effects, as shown next. However, a 
word of caution is necessary. Since the simple addition 
of the partial conductances does not reproduce the total 
conductance when there is overlap between the molecu- 
lar levels, we will not discuss the details of the Gi's, as 
they do not, by themselves, describe an experimentally 
observable quantity. Obviously, when there is no over- 
lap, as is the case for ti and ^4 = (see section IV), 
the partial conductance of each orbital is identical to the 
total conductance. 

Figure[4]shows, in the main panel, LDECA partial con- 
ductances Gi (short-dashed (green) line) and G2 (thick 
solid (red) line), for molecular orbitals | Vi) and | ^/j2), 
respectively, the total conductance Gt (long-dashed (ma- 
genta) line), which takes in account all 3 orbitals, and 
G12 (thin solid (blue) hue), as obtained through Eq. (^5]) . 
where just orbitals | ipi) and | '02) are taken in account. 
The reason why G12 and Gt are so similar is because 
orbital | V'a) is at a considerably higher energy in rela- 
tion to the degenerate orbitals and therefore its contribu- 
tion to the conductance for gate potential values around 
Vg ~ 0.5U is minimal. The phase difference A(/>i2 (in 
units of tt), as a function of gate potential, is shown in 
the inset. Note that the dip in the total conductance is 
related to a value where both partial conductances 
have the same value and A(/)i2 = it (see inset). '^^ 

The conductance features described in Fig. [5] are quite 
different from the results reported by the authors in 
Ref. 112. Although their system is the same as ours (and 
the parameters are the same) , the authors of Ref. [l3 con- 
sider, for simplicity, a regime where the total molecular 
region is described by a single level impurity, which al- 
lowed them to reduce the number of bosons in the FUS- 
BMF approximation. However, as our results show, the 
details of the internal structure of the molecule are essen- 
tial to determine its transport properties. Notice that as 
t4 increases, not only the peaks shift their positions but, 
as previously mentioned, also the structure of the peaks 
changes. 

Panels (c) and (d) in Fig. [2] display slight quantitative 
discrepancies between LDECA and FUSBMF results, al- 
though the overall qualitative agreement between the two 
techniques is quite good. These discrepancies stem from 
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FIG. 4: (color online) Detail of the conductance (calculated 
with LDECA - 11 sites) around the degenerate states | ipi) 
and I V'2): Gi (short dash (green) curve) and G2 (thick solid 
(red) curve), as defined in Eq. (22), G12 (thin solid (blue) 
curve), from Eq. (|23|l . and the total conductance (involving all 
molecular orbitals) Gt (long dash (magenta) curve) from Eq. 
(9). The inset shows the phase difi^erence between electrons 
propagating through molecular orbitals | and | V'a)- See 
text for details. 



finite-size effects in the LDECA results4S LDECA calcu- 
lations for increasingly larger exactly diagonalized clus- 
ters (not shown) indicate that the LDECA results grad- 
ually approach those from FUSBMF. This convergence 
becomes slower as Tk decreases, but the LDECA and 
FUSBMF qualitatively agree for all regimes we checked. 

Note that in the limit of strong coupling between dots 
A and C {t4 ^ t^, and for U > 2t4), a two-stage Kondo 
regime (TSK) should be expected (at half- filling) f^^ii^^ In 
this regime, dot B is weakly coupled to the band through 
the Kondo resonances of quantum dots A and C, produc- 
ing a second Kondo stage, with an exponentially smaller 
characteristic energy Ttsk ^ Tk- This special regime 
will be analyzed in a future work. 



D. S=l Kondo eflFect {t4 = ta) 

In addition to interference, the degeneracy (caused by 
symmetry) has an additional effect: it causes the two de- 
generate orbitals (when occupied by one electron each) 
to develop a ferromagnetic correlation.'*^ Indeed, when 
the structure reaches the equilateral triangle symmetry 
(i.e., i4 = t3 = 0.5), the two degenerate molecular or- 
bitals are charged simultaneously. In this case, due to 
the Kondo correlation, the first two electrons enter in 
the system with parallel spins, and the system presents 
an S" = 1 Kondo effect. This can be quantitatively ap- 
preciated by calculating the total spin for the three QDs 
as a function of gate potential, as well as the individual 
occupancy of each of the three molecular orbitals. This 
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FIG. 5: (color online) LDECA results (11 sites) for conduc- 
tance (long dash (magenta) curve), same parameters as in 
Fig- SI Charge per spin for each orbital [| ^pi), thin solid (pur- 
ple), I '02), dot-dash (black), and | tp^), short dash (blue)], and 
total spin [thick solid (green)], as a function of gate potential. 
See text for details. 



is shown by the LDECA resuhs in Fig.O where, together 
with the total conductance (long dash (magenta) curve), 
the charge occupancy for each of the molecular orbitals 
is shown (thin solid (purple) curve for orbital \ tpi): dot- 
dashed (black) curve for orbital | '02), and short dash 
(blue) curve for orbital | ips))^ and the total spin in the 
three QDs (thick solid (green) curve). Notice that the 
occupancy dependence with Vg for the two degenerate 
orbitals is not identical because they couple differently 
to the leads (orbital | 1P2) couples more strongly than 
I ipi)). As mentioned above, the maximum in the value of 
the total spin {St ~ 0.7, see thick solid (green) curve) oc- 
curs when there is approximately one electron^^ in each of 
the degenerate orbitals, which couple through an effective 
ferromagnetic interactional (Note that a value of 5 = 1 
will not be obtained for such large ratios of hopping over 
Coulomb repulsion) . This spin configuration reduces the 
ground state energy by Kondo correlating the total 5 = 1 
spin with the conduction electrons. In this region of gate 
potential, the system is in the Kondo regime, which pro- 
vides a way for the electrons at the Fermi level to cross 
from QD A to QD C. However, having two interfering 
channels at their disposal, constructed from the two de- 
generate orbitals, the conductance (long dash (magenta) 
curve) possesses a very clear Fano-like antiresonance. For 
lower gate potential values {Vg « —1.75), only orbital 
I '03) is involved in electron transport and therefore the 
conductance has the usual Lorentzian shape, with maxi- 
mum value Go, and St ~ 0.4. 



E. Deeper into Kondo and molecular regimes 

As mentioned above, the parameters in this section 
were chosen to match those in Ref. [13. As expected, and 
clearly demonstrated by the LDOS's in Fig. [3 the TQD 
system for these parameters seems to be closer to the 
intermediate valence regime than to the Kondo regime. 
Also, based on the fact that ti « t^, one may question 
if the molecular orbitals are really the most appropriate 
description of the single electron properties of the sys- 
tem. In view of that, in section IV, where the effect of 
introducing a third lead will be analyzed, the parameters 
will be changed so that the system will be deeper into 
the Kondo regime (with a larger U /F than in the current 
section). To accomplish that, we will choose U — 0.5 and 
ti = 0.2. In addition, in the next section, we will choose 
^3 — 0.4 (with < ^4 < ^3), which brings the system more 
effectively into the molecular regime (as t^/ti = 2). To il- 
lustrate both points, Fig.[S]shows the same LDOS results 
as in Fig. [3l but now for the new parameter set. It is ap- 
parent that the Kondo peaks for the new parameters are 
more well defined. For example, compare the dark solid 
(black) curves in panel (b) of both figures, which display 
the Kondo peak for the non-bonding orbital (| 1^2))- The 
peak in Fig. [H^b) clearly shows a sharper structure at the 
Fermi energy than the one in Fig. [3tb), indicating that 
this system is deeper into the Kondo regime. It is also 
apparent that the LDOS of the different molecular or- 
bitals have much less overlap in Fig. [SI underscoring the 
fact that, for the parameters to be used in section IV, the 
molecular orbitals provide a more suitable description of 
the TQD system. Nonetheless, notice that the molecular 
orbitals provide an appropriate framework to understand 
the results presented in section III as well, as it is clear 
that Figs. [3] and ini are qualitatively similar. 



IV. LOSS OF AMPLITUDE THROUGH A 
THIRD LEAD 

In this section, as just mentioned, we use different pa- 
rameters (^1 = 0.2, ts = 0.4, U = 0.5) from the ones used 
in section III. The objective is to have a larger value of 
U/T, and therefore move deeper into the Kondo regime 
and away from the intermediate valence. 

Based on a comparison of the results in Fig. [5] with 
those in Figs. [7] and [H a clear picture emerges of the ef- 
fect of a third lead connected to QD B [see Fig. [IJa)]. 
Using the labels defined in Fig. 1(a) for the QDs, let us 
qualitatively describe how the coherent propagation of 
electrons is affected by the additional lead. Assume that 
an electron is traveling from the left into QD A. After ar- 
riving at QD A, the electronic wave splits into two: one 
travels via QD C, and the other via QD B. The latter 
portion, on reaching QD B, will be split into two again: 
one travels away through the upper lead, while the other 
travels via QD C. We can view this process of 'electron 
loss' through lead P (the 'third' lead) as being a pro- 
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FIG. 6: (color online) Same LDECA calculations (9 sites) as 
in Fig. El but now for U = 0.5, ti = 0.2, t2 = 0.0, ts = 0.4, 
and ti — 0. Note that the Kondo peaks are better defined 
than the ones in Fig. [S] as there is less overlap of the LDOS 
from different molecular orbitals, indicating that for these pa- 
rameters the system is deeper into the Kondo and molecular 
regimes. Nonetheless, the qualitative similarities with Fig. [3] 
are evident, indicating that the molecular orbitals are appro- 
priate for the description of the results in section III. 
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cess of 'amplitude leakage', like that occurring at a beam 
splitter. The remaining two travehng waves (travehng 
through the triangle, in the direction of QD C) are coher- 
ent and will interfere when they propagate out of the sys- 
tem through the right lead. It will be shown below that 
the introduction of lead P does not make the electron 
propagation incoherent, since the propagation through 
overlapping molecular orbital levels clearly shows signs of 
interference, the same way as observed for t2 — 0, when 
the third lead is absent, as was discussed in Figs. [2] and 
m To analyze the results for conductance and LDOS, we 
use again the molecular orbital basis. The strategy for 
this analysis can be summarized by the following two ob- 
servations: First, by analyzing the conductance through 
each molecular orbital, one realizes that the percentage 
of the traveling wave lost through lead P will depend on 
the couphng of each molecular orbital to it. This 'loss' 
through lead P will result in a lower partial conductance 
through the molecular orbital in question. Note that, for 
a fixed value of t2 , the coupling to lead P depends only on 
the coefficient of QD B in each molecular orbital, which 
varies with the ratio ^4/^3. Second, if we assume that 
the transport through each molecular orbital is coher- 
ent (even after coupling QD B to lead P), the transport 
through two overlapping molecular orbital levels should 
give origin to interference effects, as in the case where 
lead P is not present (see section III). We will show evi- 
dence below that this is indeed the case. 



FIG. 7: (color online) Conductance as function of gate po- 
tential, Vg = VgA — VgB = Vgc , ioT different values of t2, 
showing the effect of the 'amplitude loss' due to the presence 
of the third lead. The other parameters are U = 0.5, ti = 0.2, 
ta = 0.4, and t4 = 0. Panels (a) and (b) refer to FUSBMF 
and ECA (6 sites) results, respectively. 



A. TQD in series {t4 = 0) 

Before presenting the results, we should point out that 
most of the embedded cluster results in this section were 
obtained with ECA, not LDECA.^^ As wiU be seen be- 
low, in contrast to Fig. [^^d), where the LDECA results 
suffer from minor finite-size effects, no such effects were 
detected after the third lead is connected. The higher 
symmetry obtained when t2 — ti leads to no discernible 
finite-size effects for — t^. There are two main reasons 
for that. First, the third contact provides a way for the 
'frozen' spin (present for t2 = O)**^ to delocalize from QD 
B. Second, when = and t2 = ti, any net spin would 
be equally distributed among the three QDs, diminish- 
ing its ability to suppress the Kondo effect in an ECA 
calculation.'-^ 

Let us start by turning on the connection of QD B to 
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FIG. 8: (color online) (a) to (d): Conductance as function 
of gate potential, Vg — VgA = Vgs = VgC, obtained with 
FUSBMF ['+' signs (red)] and ECA (6 sites) [stars (blue)] for 
several values of ti. Panels (a), (b), (c) and (d) show how 
the peak structure evolves from three to two peaks as the 
system goes from isosceles to equilateral triangle symmetry 
\ti = = 0.4, panel (d)]. For all panels, U — 0.5, t2 — ti — 
0.2, and ts = 0.4. Panels (e) to (g) show ECA (6 sites) results 
for the partial conductances Gi (Gi, thick solid (red) curve; 
G2, solid (blue) dots curve; G3, open squares (green) curve), 
together with the total conductance Gt (thin solid (black) 
curve). Panel (h) shows the evolution with t4 of the phase 
difference A(/!>i2 between molecular orbitals ] tpi) and \ ^2), m 
units of TT. 



lead P, by varying t2 from zero to ti = 0.2. To facilitate 
the analysis, we start with i4 = (three QDs in series, 
see Fig. [7]). In this case, the molecular orbitals are: 



^i) - [1,-a/2,1]/2 
^2) = [-1,0,1]/V2 
^3) = [l,\/2,l]/2. 



(25a) 
(25b) 
(25c) 



with El = Vg- V2t3, E2 = Vg, and E3 ^ Vg + V^ts. 
Since state | ■02) does not involve QD B, the processes 



of wave-splitting and loss of amplitude of the propagat- 
ing wave through lead P will not occur when | '1P2) is the 
state near the Fermi energy (i.e., Vg — —11/2). This re- 
sults in the partial conductance G2 through level | V'2) 
being unitary, i.e., G2 = Gq = 2e^/ft., for any value of 
t2. This is clearly what happens to the central peak in 
Fig. [71 which is associated to the molecular orbital | ?/'2), 
as previously discussed in Figs. HJa) and[31[b). The con- 
ductance through the other two molecular levels (right- 
most and leftmost peaks in Fig. [71), as mentioned above, 
will depend on the weight of QD B in | ipi) and \ ips)- 
Note that, as mentioned above, because of the choice of 
parameters {t^/ti = 2), the total conductance is basically 
the direct sum of the partial conductances when t4 — 0.0 
(as there is minimal overlap between the molecular or- 
bitals). For t4 = (see Eqs. 25 above), QD B has the 
same coefficient in | tfji) and | -03), therefore Gi = G3 for 
any value of t2 (see the identical leftmost and rightmost 
peaks in Fig. [TJ. The simultaneous and drastic decrease 
of Gi and G3 as t2 increases comes from the increase 
of the coupling of QD B to lead P, which increases the 
amplitude loss through the third lead. Note the agree- 
ment between FUSBMF (top panel in Fig. [71), and ECA 
(bottom panel) 



B. Finite t4, 

An interesting picture emerges for finite 14. In Fig. [51 
the top 4 panels show a comparison of conductance re- 
sults calculated with FUSBMF (red) signs) and ECA 
[stars (blue)] for = 0.1, 0.2, 0.3, and 0.4. The most 
salient feature in these results is the abrupt suppression 
of conductance when varies from 0.3 to 0.4. An expla- 
nation of this abrupt suppression is presented in three of 
the lower panels, which show ECA results for the total 
conductance (for t4 — 0.1, 0.3, and 0.4, in panels (e), (f), 
and (g), respectively), as well as partial conductances 
through each molecular orbital. In addition, panel (h) 
shows phase difference results [see Eq. (24)] for paths go- 
ing through either molecular orbital | ipi) or | ^2), for 
varying values of t4. 

As t4 increases, the coefficient of QD B in | ipi) in- 
creases monotonically (in absolute value) , until it reaches 
— 2/\/6 for t4 ~ t3, while it decreases monotonically for 
orbital | V'3), reaching 1/^/3 for ^4 = [see Eqs. (6) and 
(7)]. In accordance to that, Gi decreases as t4 increases, 
while G3 increases, because of the associated changes in 
the coupling of states | ipi) and | tp^) to lead P: more 
couphng (I Vi)), more 'leakage'; less couphng (| ips)), 
less 'leakage'. Note that, in panels (e) to (g) in Fig. [H 
this is evident, as the solid (red) curve corresponds to 
Gi and the open squares (green) curve corresponds to 
G3. In addition, since molecular orbital | 1P2) is inde- 
pendent of t4 [see Eq. 6(b)], the maximum value of G2 
(solid (blue) dots) is Gq for all values of t4 (no coupling 
of lead P to | 02) results in no 'leakage'). This can also 
be clearly seen in panels (e) to (g) of Fig. [51 where results 
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for G2 are shown with sohd dots (blue) . Finally, as ^4 ap- 
proaches ^3, molecular orbitals | ijji) and | 1P2) approach 
each other [see Fig. IHb)], allowing interference between 
them to strongly influence the transport properties of the 
system for Vg values where these orbitals are close to the 
Fermi energy. It so happens that the phase difference 
between the paths through these two orbitals [as calcu- 
lated in accordance to section III, Eq. (24), and shown 
in Fig. HKh)], for the relevant values of Vg, changes from 
approximately zero (for ^4 = 0.1, dashed (red) curve) to 
TT (for <4 = is = 0.4, double-dot-dash (black) curve). As 
discussed in section III, the interference will have notice- 
able effects only when the molecular orbitals | ipi) and 
I V2) are close enough in energy (for fa t^). At this 
point, the interference will be mostly destructive, as the 
phase difference is ~ tt (the resulting total conductance 
Gt is shown for all values of t4 in panels (e) to (g) as a 
thin solid (black) curve). We see the abrupt suppression 
of the central peak (associated to G2), due to its inter- 
ference with the right-side peak (associated to Gi) for 
^4 = ^3 (see panel (g) in Fig. [5]). 



One may ask why the suppression of the conductance 
in Fig. ^d) is less severe than the one in Fig. [2Kd) [note 
that there are no Fano anti-resonances in Fig. Efd)] . The 
reason is that, because of the presence of lead P, Gi is 
considerably less than the unitary conductance value Go 
[see thick solid (red) curve in Fig. Eljg)]. Therefore, de- 
structive interference cannot be total (even if the phase 
difference is tt), as G2 and Gi have widely different val- 
ues. This is not the case in Fig. [H where both partial 
conductances have similar values (being exactly the same 
at one Vg value), as in that case lead P is not present. 
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FIG. 9; (color online) LDOS calculated using LDECA (6 
sites) for U = 0.5, ti = 0.2, = 0.4, t4 — 0, and varying 
values of t2, at the particle-hole symmetric point {Vg/U = 
—0.5). Panel (a) contains results for QDs A and C and panel 
(b) for QD B. In panel (a), it is evident that coupling QD B 
to lead P {t2 = 0.1, dashed (red) curve, and t2 = 0.2, solid 
(blue) dots curve) has no effect over the density of states of 
QDs A and C, as curves are indistinguishable from the t2 = 
case (solid (black) curve). Thus, the conductance does not 
change with t2 (see Fig. [7|. Note in panel (b) that QD B 
has no spectral weight at a; = for any value of t2, implying 
that when charge transport occurs through orbital j 1P2), the 
conductance of the TQD is not affected by its coupling to lead 



CONCLUSIONS 



The results just described for the conductance of the 
central peak (| V'2)) in Fig. [3 where ^2 takes values 0, 0.1, 
and 0.2, can also be understood in terms of the LDECA 
density of states. Figure [H] shows the LDOS for QDs A 
and C [panel (a)], and QD B [panel (b)], for Vg = -C//2 
(corresponding to the central peak in Fig. [7]). Notice that 
there is no sizable change in the value of the density of 
states at the Fermi energy (w = 0), for any of the QDs, 
as ^2 varies. Since in this case the conductance is directly 
proportional to the density of states at the Fermi energy, 
this leads to a central peak in Fig.[7]that does not change 
with t2- Obviously, this independence from t2 comes from 
the fact that the density of states of QD B is very small 
in a broad interval around the Fermi energy when the 
charge transport occurs through orbital | ■02), resulting 
in lead P being effectively disconnected from the TQD 
for this value of Vg. Notice that as lead P couples to 
the other molecular orbitals (| i>i) and | V^s)); there is 
no longer a simple proportionality relation between the 
LDOS and the conductance^^ so that no simple direct 
connection can be made between the LDOS at the Fermi 
energy and the conductance (as was done in Fig. [31). 



We have studied the transport properties of a TQD 
in the molecular regime coupled to leads. By applying 
the FUSBMF, and the EGA and LDEGA approaches, 
we have calculated the conductance of the TQD system 
for different symmetries and different configurations of 
the leads. For the two-leads case, we have calculated the 
conductance for both series and triangle configurations. 
In the series configuration, LDEGA and FUSBMF results 
agree with each other and with the results for the molecu- 
lar regime obtained in Ref.[l^, where the Kondo effect has 
been studied in detail. In the triangular symmetry, the 
quantitative results obtained by FUSBMF and LDEGA 
differ slightly as the equilateral symmetry is approached, 
due to a finite-size effect in the relatively small clusters 
accessible to LDEGA, although agreement is still very 
good [see Fig. [2Kd)]. The suppression of conductance in 
the regime where approximately two electrons occupy the 
triangle was explained by LDEGA as an interference ef- 
fect between two degenerate molecular orbitals, utilizing 
the concept of partial conductance. In addition, our re- 
sults for triangular symmetry differ from those presented 
recently in Ref. 17. We believe that the approach pur- 
sued here, where details of the internal structure of the 
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interacting region of the system are taken fully into ac- 
count, are very important to explain the conductance of 
the TQD system. In fact, our results show that changes 
in the internal couplings of the TQD dramatically change 
the features of the conductance. We also found that the 
degeneracy of the molecular orbitals at equilateral sym- 
metry, when two electrons occupy the TQD, induces an 
effective ferromagnetic interaction between the spins lo- 
calized in the interacting region;^ leading to an 5* = 1 
Kondo effect. 

In the TQD series configuration, our results show that 
the third lead produces a strong suppression in the bond- 
ing and anti-bonding orbital conductance peaks (Fig. [7]). 
The non-bonding peak, however, remains unchanged, 
since this orbital does not have the appropriate symme- 
try to couple to lead P. This suppression of conductance 
can be seen as a 'loss of amplitude' through lead P, sim- 
ilar to the effect occurring with beam splitters in optics. 
If one thinks of the conductance in terms of transmis- 
sion of waves through the interacting region, the intro- 
duction of lead P provides an additional transmission 
channel, which clearly affects the conductance between 
leads L and R. This 'loss of amplitude' idea is then used 
to understand the conductance results in the triangular 



symmetry. In particular, it explains why the interfer- 
ence effects seem less effective in suppressing the con- 
ductance in the equilateral symmetry (i.e., why no Fano 
anti- resonance occurs): the 'loss of amplitude' prevents 
the conductance through molecular orbital | ipi) from 
reaching the unitary limit, leading to a decrease in the 
destructive interference, as discussed in detail in Fig. [8l 
We should remark that the excellent overall quantitative 
agreement of results obtained with FUSBMF, EGA, and 
LDECA"*-"^ (which rely on totally different approxima- 
tions) makes our conclusions much more reliable and ro- 
bust. Moreover, the combination of techniques allows a 
better insight into the physics of the different geometries. 
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